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Prepared  by: 

R.  E.  Wilson 

ABSTRACT:  A  momentum-integral  method  for  calculating  the  boundary- 
layer  growth  on  slightly -blunted  cones  has  been  given  in  a  previous  paper 
by  the  writer.  The  method,  which  was  applied  to  cones  with  fully 
laminar  flow,  is  extended  in  the  present  paper  to  turbulent  flow.  In 
addition:  (1)  transition  data  on  cones  at  hypersonic  speeds  are  examined 
to  show  the  great  eifect  of  slight  nose  blunting  on  interpretation  of  the 
results,  (2)  agreement  between  measured  and  calculated  cone  drag 
coefficients  for  fully  laminar  flow  is  demonstrated,  and  (3)  turbulent 
friction  and  boundary -layer  thickness  calculations  are  shown  to  agree 
with  the  limited  amount  of  ballistics  range  data  which  are  available. 
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*  ^(pu/piUx)[l  -  (u/ux)]dy  = 


boundary -layer  momentum  thicknes  ; 


*  absolute  viscosity 

-  Reynolds  number:  Rxc  =  pc ucx/pc  ;  R0C  =  Pcuce/Mc  J 
Rbc  =  *  ^Qi  =  Pi  ui  »  ^L®  =  Pcd^osL/pob 


=  density 
=  cone  half  angle 
=  shear  stress 

=  shock  wave  angle  fcr  sharp  cone 
=  local  shock  wave  angle  for  blunt  cone 
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w  =  conditions  at  the  wall 

1  =  local  conditions  at  the  outer  edge  of  the  boundary  layer 

•  =  freestream  conditions 

T  =  stagnation  conditions 
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LAMINAR  AND  TURBULENT  BOUNDARY  LAYERS 
ON  SLIGHTLY -BLUNTED  CONES  AT 
HYPERSONIC  SPEEDS 


INTRODUCTION 


The  momentum -Integral  method  presented  in  reference  (1)  yields 
local  flow  conditions  over  slightly -blunted  cones  and  at  the  same  time 
y.eids  values  for  the  boundary-layer  thicknesses,  6,  6*, and  6.  In 
addition,  shear  stress  and  heat-transfer  distributions  over  the  cone  are 
determined.  Reference  (1)  gives  a  skin-friction  relation  and  expressions 
for  the  boundary -layer  thicknesses  which  apply  only  to  laminar  flow.  In 
the  present  paper  relations  are  given  for  turbulent  boundary-layer  flow  so 
that  the  solution  of  the  momentum-integral  equation  can  be  continued  after 
transition. 

For  the  laminar  boundary -layer  case,  references  (2)  and  (3)  give  a 
method  for  calculating  the  local  Mach  number  variation  along  the  cone 
which  assumes  local  similarity  and  makes  use  of  a  mass  flow  function 
of  reference  (4).  The  method  does  not  account  for  the  effect  of  Mach 
number  gradient  on  boundary-layer  growth.  In  addition,  the  method 
holds  only  for  the  case  of  highly -cooled  walls.  The  momentum -intec,  il 
method  presented  in  reference  (1)  and  used  here  does  not  have  these 
limitations.  A  comparison  was  made  in  reference  (1)  between  the 
momentum -integral  method  and  the  method  of  reference  (2) .  There  is  a 
difference  between  references  (2)  and  (3)  which  greatly  affects  numerical 
results  and  therefore  both  references  have  been  compared  here  'vith  the 
momentum -integral  method. 

The  shear  stress  at  the  outer  edge  of  the  boundary  layer  is  neglected 
in  the  momentum-integral  equation.  The  accuracy  in  this  assumption  is 
commented  on  in  reference  (1).  An  expression  for  calculating  the  ratio, 
Ti/tw  for  the  laminar  boundary  layer  is  derived  here  so  that  this  assump¬ 
tion  can  be  checked  for  any  specific  case. 

The  ballistics  ranges  at  the  U.  S.  Naval  Ordnance  Laboratory  have 
been  used  to  measure  transition  on  cones,  the  drag  of  fully  laminar 
cones,  the  drag  of  cones  with  turbulent  flow,  and  even  to  get  some  idea  of 
the  thickness  of  the  turbulent  boundary  layer  at  the  base  of  a  cone .  These 
data  have  been  examined  and  pertinent  calculations  using  the  momentum- 
integral  method  have  been  made.  The  calculations  have  been  compared 
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with  experiment  in  some  cases  and  in  others  have  been  used  in 
analyzing  the  experimental  data. 


Theoretical  Results 


Using  the  coordinate  system  on  figure  1,  the  momentum -integral 
equation  for  the  conical  part  of  a  slightly -blunted  cone  is  given  in 
reference  (1)  as 


The  assumptions  made  in  deriving  equation  (1)  are:  (1)  a  constant 
pressure  on  the  conical  surface  equal  to  the  inviscid  sharp  cone  pressure, 
(2)  a  thermally  and  calorifically  perfect  gas,  (3)  adiabatic  flow  outside 
the  boundary  layer,  and  (4)  2t1/piU?  «  Cf. 


Boundary-Layer  Parameters 


The  expressions  for  Cf,  6/0,  and  H  are  taken  from  the  flat  plate 
case  and  evaluated  at  the  local  flow  conditions  at  each  point  along  the 
body.  The  expressions  which  apply  when  the  boundary  layer  is  iaminar 
are  given  in  reference  (!).  The  expressions  which  will  be  used  when 
the  boundary  layer  is  turbulent  are  taken  from  the  flat  plate  results 
given  in  reference  (5)  and  discussed  briefly  below. 


The  flat  plate  analysis  is  based  on  von  K£rm2n's  similarity  theory 
for  incompressible  flow.  An  apparent  Prandtl  number  near  unity  is 
assumed  so  that  the  density  can  be  approximated  by 


where 


Using  equation  (2)  and  taking  a  constant  wall  temperature  along  the  plate, 
a  derivat.on  similar  to  the  insulated  plate  derivation  of  reference  (6) 
was  carr  i.ed  out  to  obtain  the  skin  friction  relations. 


2 
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(sln“x  r,  -  sin"1  rjf)  OcfTTi/T1)‘^=  -6.43  +  logw|~£- Jfl  CfR^| 

0.242  (sin*1  rj  -  sirf1  rjp)  (/3cfTTi/t  1)‘^=  -1.968  +  loglc^^  —  cfrJ 

Taking  z  here  as  the  distance  from  the  plate  leading  edge,  D  as  the 
plate  drag,  and  as  the  velocity  at  the  edye  of  the  laminar  sublayer 

R  =  pi  ul  z/ Hi 

Cf  =  2  Tw/piUia  «  local  skin  friction  coefficient 
cF  =  2D/piUxaz  =  mean  skin  friction  coefficient 

s  =  u*(Tw/pw)* 


(3) 


(4) 


Experimental  values  of  the  sublayer  parameter,  S,  have  been  collected 
In  reference  (7)  where  it  Is  shown  that  a  good  fit  to  the  data  is  given  by 

S  =11.5  + 6.6  (Tei  -  Tw)/Tei  (5) 


The  recovery  temperature,  Tel ,  for  the  turbulent  boundary  layer  is 
given  by 


=  (1  +  Pr* 


M=  )T, 


(6) 


The  functions  t),  Tjp  and  rjF  are  defined  by 


n  =|20  -  (ttx  -  Tw)/TTi }/  $ 

«  (If-  $  -  <tt,  -  tw)/Tt,}/* 

if  f]-  frr,  -  WTtJ/* 
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where 


*  .{4  e[l  -  (TXi  -  Tw)/TtJ+  [(TT>  -  Tw)/TTl]a}i 

Theoretical  results  based  on  von  K£rmcLn‘s  similarity  theory  are 
given  in  both  references  (7)  and  (8) .  The  relations  given  in  both  references 
differ  somewhat  from  equations  (3)  and  (4).  The  present  results  have  been 
compared  with  experimental  skin  friction  coefficients  of  references  (6), 

(9),  (10),  (11),  (12),  and  (13)  which  cover  the  range  of  Mx  from  1.7  to  7.0 
and  (TTj  -  Tw)/Ttx  from  0.04  to  0.84.  There  is  reasonable  agreement 
between  theory  and  experiment.  It  should  be  noted  that  equations  (3)  and 
(4)  and  the  expressions  for  17,  rjf ,  tip,  and  differ  slightly  from  those 
given  in  AGARDograph  97.  Where  the  total  temperature  Tp  appears 
here,  the  recovery  temperature  Tei  appears  in  AGARDograph  97.  The 
change  was  made  after  comparing  the  analytical  results  with  high  Mach 
number  friction  drag  data  recently  obtained  at  the  Naval  Ordnance 
Laboratory.  With  the  assumption  Pr  =  1,  it  is  in  fact  Tp  which  appears 
in  the  analysis .  The  previous  use  of  Te  was  an  attempt  £0  account  for 
the  fact  that  Pr  ^  1 . 


To  integrate  equation  (1)  over  a  cone,  it  is  necessary  to  obtain  cf 
as  a  function  of  Rq  .  It  is  assumed  that  the  expression  derived  below 
for  a  flat  plate  applies  here .  For  a  flat  plate 


cF  =  2  0/z  =2  Rei/R 
Cf=  d(cFR)/dR 


With  equations  (4),  (7),  and  (8)  it  can  be  shown  that 


cf  =  Cp 


-1.968  + log 


101 


El. 

Mw 


R 


e.)  +°-2«  {z^M-4) 


-1.099  + log 


><$r  ^Rei)  +  0-242  (2  3  wJ/4'4) 


With  equation  (7),  equation  (4)  can  be  written 


(7) 

(8) 


(9) 


0.242  (sin-1  rj  -  sin  -1 


-1.968 


El 

Mw 


In  effect,  cF  can  now  be  eliminated  between  equations  (9)  and  (10). 
Since  this  can  only  be  done  numerically,  equatio;.  (9)  and  (10)  serve  as 


4 


"  MM 
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the  relation  between  c{  and  R0i  •  The  coefficient  cp  is  merely  a 
parameter  which  is  eliminated  numerically  when  making  computations. 

The  parameters  6/0  and  H  =  6*/0  have  been  approximated  by 
assuming  that  the  boundary -layer  velocity  distribution  is  given  by 


*  -W 


Values  of  6/6  and  H  have  been  calculated  using  equations  (2)  and  (11)  and 
the  definitions  for  6*  and  0.  The  results  can  be  approximated  by 


6  _  (n  + 


6  n  +  2  (Tg^  JIwj/Tg^ 


+  y  ~1ul  {l -34  +  1.38 


TTx  - 


n  +  2 


TTx-T 

TTi 


M?  |l  4-  1.20  (l  -  Tl^  Tw^  j.  (13) 


Values  of  6/6  and  H  are  not  very  sensitive  to  n  which,  from  experiment, 
varies  slowly  with  Reynolds  number.  For  a  nominal  value  of  n  =  7,  a 
range  of  Mx  from  0  to  95  and  ^Tg^  -Tw)/  Ttx  from  0  to  1,  the  agreement 

between  the  numerical  results  and  equations  (12)  and  (13)  is  generally 
better  than  3%  and  1.5%,  respectively. 

The  derivative  d  £  (6/6)  -  H  ]  /dMx  will  be  required  later.  For 
Tw  =  constant  and  n  =  constant,  the  following  result  can  be  obtained  from 
equations  (12)  and  (13) 

(y-l)Mx  0.34  +  I  Tw^  (14) 

Mach  Number  Gradient 

Reference  (1)  presents  two  relations  between  the  local  Mach  number, 
Mx ,  at  the  outer  edge  of  the  boundary  layer  and  the  local  shock  wave 
angle,  The  first  relation  is  derived  by  equating  the  mass  flow  in  the 
boundary  layer  to  that  in  a  stream  tube  ahead  of  the  bow  shock.  The  bow 
shock  is  assumed  to  be  hyperbolic  and  described  by 

a3/r|  =  tan3  cos  -  tansu;  (15) 


where  the  shock  curvature  on  the  axis,  a  ,  is  related  to  the  body  nose 
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radius  by  a  M  1 . 5b.  The  relation  is 


(16) 


The  second  relation  assumes  isentroplc  flow  along  each  streamline  from 
the  bow  shock  to  the  outer  edge  of  the  boundary  layer .  With  this  assump¬ 
tion  and  the  oblique  shock  equations  for  a  perfect  gas 

y  ♦  1 _ 

2yM^sin3ws-(y-l) 


(y»l)&£slna<i;s  [2  +  (y-l)M^] 
2  [2  +  (y-l)M3  sin3 


(17) 


Making  use  of  the  assumption  that  Px  *  constant  on  the  cone,  differentiat¬ 
ing  equations  (16)  and  (17)  with  respect  to  x  and  eliminating  do>s/dx 
results  in  an  expression  which  can  be  used  for  calculating  the  Mach  number 
gradient,  dMx/dx. 


26 _  |2  +  (y-l)Mi  -  r(6/e)-Hl)dM1  =  d8  +  6  dr 

Mx  [  2  +  (y-l)Mg|  Fx  j  dx  dx  r  dx 


(18) 


where 


Mx  [2^(y-l)M?]  dC(6/e)-H] 

F,-  -  -{2.(y-l)Mj-(|  -  H)}/{l  ♦  (y-l)M?  *  2  [(8/e)  _H]  dMi 


|2yM*slnan)s-(y-l)]f2  .  (y-l)Mj,slnaais|) 
-  tnn8^  (Ua  Sin3  tAlr*  -  f 


(19) 


M3tan3oos  sec3a^ 

2  (tan3  ojs  -  tan3cu)  (M3,,,  sin3Ws  -  l)a 


For  the  laminar  boundary  layer  case  treated  in  reference  (1),  the  difference 
(6/0)  -  H  which  appears  in  equations  (16),  (18),  and  (19)  above  (and  also  in 
equation  (1))  is  constant.  For  the  turbulent  case  this  is  not  true  and  compli¬ 
cates  the  calculations.  Equations  (12)  and  (13)  must  be-used  to^alculate 
(6/0)  -  H  and  equation  (14)  must  be  used  to  calculate  d  (6/6)  -Hj  AjMx  which 
appears  in  equation  (19) .  L  -* 
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Integration  of  the  Momentum  Equation 

Making  use  of  the  definitions  for  R.  and  R  ,  equations  (1)  and  (l£) 
yield  be  xc 


Integration  of  equation  (20)  for  the  laminar  case  is  treated  in  reference  (3). 
For  the  turbulent  case,  the  integration  starts  at  the  transition  point  where 
Rxc  =  Rxctr  and  Rqc  =  Rectr*  The  integral  form  of  equation  (20)  is  then 


In  order  to  obtain  Rc„  as  a  function  of  Rv/%  from  equation  (21),  the  value 

we  xc 

of  eitner  Rgctr  or  Rxctr  must  known  or  assumed.  The  other  value  is 

obtained  by  integrating  equation  (20)  for  the  laminar  portion  of  the  boundary 
layer  using  reference  (1).  To  proceed  with  the  integration,  equation  (21) 
is  used  v/ith  Cf  given  by  equations  (9)  and  (10)  and  Fx  given  by  equation  (19) . 
It  should  be  noted  that  to  obtain  Cf  as  a  function  of  R^,  put 

Rqx  =  piU1/icRgc/pcuc^i1  in  equations  (9)  and  (10).  With  constant  pressure 
on  the  cone  and  adiabatic  flow  outside  the  boundary  layer, 

PcUc/PiUx  =  McTf/MxTj  and  Tx/Tc  *[2  +  (y-l)Mj]/[2  +  (y-l)M?  ]  • 

As  equation  (21)  is  integrated,  the  value  of  Mx  is  obtained  along  the  cone 
by  finding  the  values  of  Mx  and  u:s  at  each  point  which  satisfy  equations 
(17)  and  (16)  which  may  be  written 


1.125  Poe  Mod 

2  +  (y-l)MS 

sin  a  PXMX 

+  (y-l)M?! 

Rbc  [(6/6)-h]  (tan3ws-tanaw) 


(22) 


The  above  method  assumes  that  6  is  continuous  at  the  transition 
point.  This  results  in  discontinuities  in  the  other  thicknesses,  6*  and  6« 
In  addition,  there  is  a  disc'  ntinuity  in  Mx  corresponding  to  the  dis¬ 
continuity  in  6.  This  results  from  the  Mach  number  gradient  normal  to 
the  surface  caused  by  the  shock  curvature.  The  value  of  Mx  for  the 
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turbulent  boundary  layer  at  the  transition  point  can  be  obtained  from 
equations  (17)  and  (22).  Put  the  transition  point  values  of  Rxr>  and  R0C 
In  equation  (22)  and  calculate  the  value  of  (6/0)  -  H  required  from 
equations  (12)  and  (13) . 

As  equation  (21)  is  integrated  to  obtain  R~  and  local  properties 
along  the  cone,  the  local  skin  friction  is  determined  at  each  point  and 
all  the  thicknesses  (8,  6*,  and  6)  can  be  calculated.  As  in  reference  (1), 
the  heat-transfer  rate  can  be  calculated  assuming  that  the  Stanton 
number  is  given  by  the  flat  plate  Reynolds  analogy.  For  the  turbulent 
boundary  layer  the  local  recovery  temperature  must  be  determined  for 
each  point  from  equation  (6) . 

It  should  be  noted  that  equation  (21)  applies  to  the  case  of  a  sharp 
cone  by  putting  Fx  =  0. 


Numerical  and  Experimental  Results 

Figure  2  presents  numerical  results  showing  the  variation  of 
local  Mach  number  along  an  8°  half  angle  cone  with  three  bluntnesses . 
The  flow  passing  through  the  curved  portion  of  the  bow  shock  forms  a 
variable  entropy  layer  over  the  cone.  As  the  boundary  layer  swallows 
this  variable  entropy  layer,  the  local  Mach  number  increases  and 
approaches  the  sharp  cone  value.  The  solid  curves,  taken  from 
reference  (1),  shew  the  effect  of  nose  radius  Reynolds  number,  R^c, 
on  the  laminar  boundary  layer.  Experimental  results  discussed  later 
indicate  that  the  momentum  thickness  Reynolds  number  based  on  local 
properties,  R01 ,  may  be  approximately  800  at  transition.  Assuming 
no  dependence  on  Mach  number  and  taking  this  value,  the  dotted  curves 
on  figure  2  show  the  great  effect  of  the  turbulent  boundary  layer  on  local 
Mach  number.  For  the  lowest  value  of  R^c,  the  sharp  cone  Mach 
number  is  reached  before  transition  occurs.  The  aft  movement  of 
transition  with  increasing  R^c  is  due  to  the  decreasing  local  unit 
Reynolds  number  and  the  assumption  that  Rgi  at  transition  is  constant. 

References  (2)  and  (3)  give  a  method  for  calculating  the  local  Mach 
number  variation  for  the  laminar  boundary  layer.  Instead  of  integrating 
the  momentum  equation  as  is  done  here,  the  method  assumes  local 
boundary  layer  similarity  and  makes  use  of  a  mass  flow  function  given  by 
reference  (4).  The  momentum -integral  results  from  reference  (1)  are 
compared  on  figure  3  with  results  using  references  (2)  and  (3) .  The  bow 
shock  given  by  equation  (15)  was  used  for  all  calculations.  There  is  a 
great  difference  between  the  results  of  references  (1)  and  (2). 
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References  (1)  and  (3)  are  in  good  agreement  for  the  case  calculated. 

In  reference  (2)  the  ratio  of  local  stagnation  pressure  to  Pitot  pressure 
was  Included  in  defining  a  transformed  coordinate  along  the  body.  There 
seemed  no  justification  for  Including  this  ratio  and  it  was  omitted  when 
the  same  authors  prepared  reference  (3). 

The  effect  of  a  slight  bluntness  on  local  conditions  can  significantly 
affect  the  interpretation  of  experimental  data.  The  interpretation  of 
transition  data  from  ballistics  range  firings  is  a  pertinent  example. 
Figure  4  taken  from  reference  (14)  presents  drag  coefficient  data  for  a 
slightly-blunted  cone.  As  the  freestream  Reynolds  number  is  increased, 
transition  moves  upstream  in  the  wake.  At  some  Reynolds  number  the 
wake  will  be  completely  turbulent  and  transition  will  be  located  at  the 
base  of  the  body.  With  a  further  increase  in  Reynolds  number,  transition 
moves  forward  on  the  body  and  the  drag  rises.  The  open  symbols  on 
figure  4  indicate  that  the  base  flow  was  laminar  on  shadow  photographs. 
The  solid  symbols  indicate  that  the  base  flow  was  completely  turbulent. 

At  the  freestream  Reynolds  number  for  which  the  drag  rise  begins, 
transition  occurs  at  the  base .  The  drag  curve  can  thus  be  used  to 
determine  a  transition  Reynolds  number  for  slender  cones.  This  is  a 
useful  technique  at  hypersonic  speeds  wrhen  the  shock  is  so  close  to  the 
surface  that  the  boundary  layer  is  difficult  to  see  on  a  schlieren  or 
shadow  photograph.  When  the  cone  is  slightly  blunted,  a  difficulty  now 
arises  in  stating  the  local  Mach  number  and  Reynolds  number  at  which 
transition  occurs.  In  order  to  resolve  this  difficulty,  the  boundary -layer 
growth  over  the  body  must  be  calculated.  Laminar  flow  only  need  be 
considered  since  transition  occurs  at  the  base.  Values  of  R^c  at 
transition  of  3.8  x  10*  and  6.3  x  104  were  computed  at  Mach  numbers  of 
9  and  13.3,  respectively,  for  the  data  shown  on  figure  4.  For  these 
values  of  R^c,  the  boundary  layer  growth  over  the  cone  was  computed 
using  the  method  of  reference  (1).  The  variations  in  Mach  number  over 
the  cone  are  shown  on  figure  5.  The  arrows  marked  transition  indicated 
the  values  of  Rxc  at  the  base  of  the  cone.  It  is  evident  that  the  swallowing 
process  had  not  been  completed  at  either  Mach  number  9  or  13.3  when 
transition  occurred.  Using  the  boundary-layer  calculations,  the  actual 
transition  conditions  can  be  compared  with  those  which  would  be  obtained 
if  the  bluntness  were  neglected  and  the  cone  assumed  to  be  sharp. 

Figure  6  shows  the  transition  Reynolds  number  based  on  length  along  the 
cone  as  a  function  of  Mach  number.  In  one  case  tie  Reynolds  number  is 
based  on  local  conditions  at  the  base  and  the  local  Mach  number  is 
plotted  on  the  abscissa.  This  is  compared  with  the  base  where  Reynolds 
number  based  on  sharp  cone  conditions  is  plotted  vs.  the  sharp  cone 
Mach  number.  The  magnitudes  of  the  results  are  quite  different  and, 
in  addition,  the  trend  of  transition  Reynolds  number  with  Mach  number 


9 


NOLTR  66-54 


Is  different.  It  Is  perhaps  more  pertinent  to  plot  values  of  R«  vs.  M 
at  transition  for  the  blunt  and  sharp  cone.  These  data  are  shown  on 
figure  7  and  the  comment  made  for  figure  6  can  be  repeated.  At  a  Mach 
number  of  7,  the  data  indicate  that  the  local  Reynolds  number,  Rqj  ,  at 
transition  is  about  800. 

Slight  blunting  can  significantly  affect  the  local  shear  stress. 
Calculations  for  the  laminar  boundary  layer  taken  from  reference  (1)  are 
shown  on  figure  8.  The  ratio  Cfc/c|  is  equal  to  the  ratio  of  the  wall 
shear  stress  on  the  blunted  cone  to  that  on  a  sharp  cone,  since  both  skin 
friction  coefficients  are  based  on  the  sharp  cone  dynamic  pressure.  The 
results  of  figure  8  show  that  the  bluntness  should  be  taken  into  account 
when  calculating  the  friction  drag.  Reference  (14)  presents  experimental 
and  calculated  values  of  total  drag  coefficient  for  fully  laminar  cones. 

The  data  are  plotted  in  reference  (14)  vs.  the  hypersonic  interaction 
parameter.  They  are  replotted  here  on  figure  9  as  a  function  of  free- 
stream  Reynolds  number.  There  is  good  agreement  between  calculation 
and  experiment.  The  calculations  included  the  pressure  drag  on  the 
cone,  the  base  drag,  the  friction  drag  calculated  using  the  method  of 
reference  (1),  an  induced  pressure  correction  to  the  cone  pressure  drag, 
and  induced  pressure  and  transverse  curvature  corrections  to  the 
friction  drag.  However,  the  difference  between  the  cone  pressure  drag 
coefficient  shown  on  figure  9  and  the  total  drag  coefficient  is  largely  due 
to  the  friction  drag  calculated  using  the  method  of  reference  (1).  This 
indicates  that  the  momentum -integral  method  for  slightly-blunted  cones 
gives  reliable  laminar  friction  drag  predictions.  It  is  of  interest  to 
note  from  figure  9  that,  when  plotted  vs.  freestream  Reynolds  number, 
the  total  drag  coefficient  seems  to  be  almost  independent  of  Mach 
number . 

In  deriving  the  momentum -integral  equation,  equation  (1),  it  was 
assumed  that  the  shear  stress  at  the  outer  edge  of  the  boundary  layer 
could  be  neglected;  i.e.,  2rl/p1Ui«Cf.  Reference  (1)  stated  that  this 
assumption  had  been  checked  for  the  numerical  results  given  on 
figure  8.  This  check  can  be  made  as  follows.  The  change  in  ut  in  the  y 
direction  between  two  adjacent  stream  tubes  at  the  outer  edge  of  the 
boundary  layer  is  equal  to  the  change  in  ux  in  the  x  direction  between  the 
two  points  at  which  these  streamtubes  enter  the  boundary  layer.  In 
addition,  the  mass  flow  between  these  two  streamtubes  is  equal  to  the 
increase  in  the  mass  flow  in  the  boundary  layer.  From  these  statements 
it  follows  that 
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The  value  of  n  Is  given  by 

Ti  ■  *@r), 

With  the  above  relations  and  equations  (18)  and  (20),  it  can  be  shown  for 
the  case  of  laminar  flow  on  a  cone,  where  (6/6)  -  H  =  constant  =  6.41,  that 


=  {3.405  =fRei[1+(y.1)M?]  +  if 

The  above  result  is  Independent  of  the  bluntness  Reynolds  number,  Rbc. 

Making  use  of  reference  (1)  to  compute  c^Rfi  and  Fx,  Tx/r  v  ha:: 

been  computed  and  plotted  vs.  Mx  on  figure  10.  This  result  applies  to 
all  three  bluntnesses  on  figure  8.  It  can  be  seen  that  the  error  in 
neglecting  Ti  is  small. 


Figure  11  shows  a  shadow  photograph  from  a  ballistics  range  shot 
of  a  9°  half  angle  cone  at  a  freestream  Mach  number  of  10.3.  The 
enlargement  shows  the  flow  in  the  base  region.  The  boundary  layer  is 
almost  completely  turbulent  and  data  which  can  be  obtained  from  the 
shot  have  been  compared  with  calculations.  The  boundary -layer  growth 
was  calculated  using  reference  (1)  for  the  laminar  region  and  the  present 
results  for  the  turbulent  region.  A  value  of  n  =  7  was  used  for  the 
velocity  distribution  of  equation  (11).  Although  the  location  of  transition 
is  impossible  to  determine  accurately,  it  appeared  from  an  enlargement 
of  figure  11  to  occur  about  20%  of  the  body  length  from  the  nose.  On 
Integrating  the  momentum -integral  equation,  it  was  found  that  at 
transition  was  approximately  300-  The  computed  variation  in  Mach 
number  at  the  outer  edge  of  the  boundary  layer  is  shown  on  figure  12. 

It  can  be  seen  that  transition  occured  at  Mi  «  3.  The  value  of  R0  in  this 
case  is  much  less  than  the  value  of  800  at  Mx«  7  given  on  figure  7.  It 
should  be  noted  that  there  was  appreciable  melting  of  the  tip  of  the  model 
shown  on  figure  11.  This  was  not  the  case  for  the  data  of  figure  4.  The 
computed  variation  of  boundary -layer  thickness,  6,  is  shown  on  figure  13. 
Values  of  6  at  the  base  were  measured  at  the  top  and  bottom  of  the  body 
from  an  enlargement  of  figure  11.  Measurements  were  also  made  on  a 
repeat  shot.  The  agreement  between  calculation  and  experiment  on 
figure  13  is  remarkable  and  probably  fortuitous  considering  the  difficulty 
in  defining  the  edge  of  the  boundary  layer.  Total  drag  measurements  were 
made  for  the  two  shots.  The  measured  drag  coefficients  of  0.0854  and 
0.0839  compare  with  a  calculated  value  of  0.0929.  The  check  between 
calculation  arid  experiment  is  reasonable.  However,  no  statement  can  be 
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made  about  agreement  between  experiment  and  calculated  friction  drag. 

This  Is  due  to  uncertainties  In  base  pressure  and  model  tip  diameter  and 
the  fact  that  the  friction  drag  was  only  about  20%  of  the  total.  At  this  point 
it  is  of  interest  to  return  to  figure  4,  where  the  check  between  calculated 
and  experimental  drag  results  is  quite  good .  At  the  higher  Reynolds  numbers 
when  considerable  turbulent  flow  exists,  the  friction  drag  is  30%  or  more 
of  the  total.  This  increase  in  the  relative  magnitude  of  the  friction  drag  is 
due  primarily  to  the  small  cone  half  angle  of  6.3°  .  The  friction  drag  when 
turbulent  flow  was  present  was  calculated  for  figure  4  assuming  a  constant 
Rqx  at  transition. 

One  of  the  interesting  features  of  figure  11  is  the  very  wide  turbulent, 
wake  in  the  base  region,  with  the  trailing  shocks  originating  deep  inside 
the  wake.  From  the  calculations  of  boundary -layer  thickness  and  also  from 
figure  11,  it  is  apparent  that  the  boundary  layer  occupies  about  half  of  the 
region  between  the  model  surface  and  the  bow  shock.  Thus  when  tho  flow 
expands  around  the  corner  at  t.he  base,  a  great  part  of  the  base  region  is 
filled  with  the  boundary -layer  flow. 
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FIG.  1  COORDINATE  SYSTEM 


FIG.  2 


EFFECT  OF  BLUNTNESS  ON  LOCAL  MACH  NUMBER 
Mo,  =  14.9,  Mc  =  10.0,  0=8°,  Too  =  Tw  =  530°R 
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FIG.  4 

COMPARISON  OF  CALCULATED  AND  EXPERIMENTAL  CONE  DRAG 
COEFFICIENTS:  0  =  6.3°,  b/B  =  0.035,  Too=Tw=540°R 


